Draft version February 2, 2008 

Preprint typeset using 1^1^^ style cmulatcapj v. 10/09/06 



O 

o 

> 

o 
O 

(N 

6 



> 
O 

o 
d 

o 



X 



INSTABILITIES OF SPIRAL SHOCKS - II. A QUASLSTEADY STATE IN THE MULTLPHASE 

INHOMOGENEOUS ISM 

Keiichi Wada^ 

National Astronomical Observatory of Japan, Mitaka, Tokyo 181-8588, Japan 
E-mail: wada.keiichi@nao.ac.jp 

Draft version February 2, 2008 

ABSTRACT 

The "galactic shocks" ()Fuiimotalll968l : [Roberts! 1 1969f ) is investigated using a full three-dimensional 
hydrodynamic simulations, taking into account self-gravity of the ISM, radiative cooling, and star 
formation followed by energy feedback from supernovae. This is an essential progress from the previous 
numerical models, in which 2-D isothermal, non-self-gravitating gas is assumed. We find that the 
classic galactic shocks appears is unstable and transient, and it shifts to a globally quasi-steady, 
inhomogeneous pattern due to non-linear development of instabilities in the disk. The spiral patterns 
consists of many GMC-like dense condensations, but those local structures are not steady, and they 
evolves into irregular spurs in the inter- arm regions. Energy feedback from supernovae do not destroy 
the quasi-steady spiral arms, but it mainly contributes to vertical motion and structures of the ISM. 
The results and methods presented here are a starting point for more consistent treatment of the 
ISM in spiral galaxies, in which effects of magnetic fields, radiative transfer, chemistry, and dynamical 
evolution of a stellar disk are taken into account. 

Subject headings: ISM: structure, kinematics and dynamics — method: numerical 



1. INTRODUCTION 

A standing spiral shock solution was first discovered in 
the 60s numerically in a rotating gas disk in the galactic 
potential with a tightly wound spiral perturbation (i.e. 
the "pitch angle" is very small). Since then in the most of 
theoretical studies on the galactic "sho ck" (e.g. Fuiimot3 
1968; Roberts igeg": 'Woodward' '197^ [Johns, fc NelswH 
1986. : .Lubow. Cowie. fc Balbusi il986l ). the interstellar 



medium (ISM) is treated as a isothermal and homoge- 
neous fluid with two-dimensional approximation, or it 
is modeled as a iV— body system of small cloudlets (e.g. 
lTomisalffllll986[ ). Gl obal evolution of the sp iral shock was 
studied in the 80s (jJohns. fc NelsonI [T98& ) using time- 
dependent, two-dimensional (2D) hydrodynamic simula- 
tions show that spiral shocks are stable and long-standing 
for various pitch angles. 

However, the steady, smoothed galactic shock does 
not consistently explain the complicated distribution of 
the ISM around spiral arms and the inter-arm substruc- 
tures akin to so-called 'spurs' or 'feather s'^ observed in 
real spiral galaxies ( Elmegrecn 1980; Sco ville et al.l[200ll : 
ICalzetti et al] 120051: ILa Vigne et all 120061 1 Moreover, 
observed molecular clouds in galactic disks do not match 
the picture of hydrodynamic shocks in a uniform media. 
It was however shown by full 2D global simulations of a 
non-self-gravitating, isothermal gas disk that the spurs 

^ Department of Astronomical Science, The Graduate Univer- 
sity for Advanced Studies, Osawa 2-21-1, Mitaka, Tokyo 181-8588, 
Japan. 

^ In this paper, we refer the terms 'spurs' or 'feathers' as inter- 
arm gas structures, which often show quasi-periodic fea tures as- 
sociat ed with main spiral ar ms. Se e Paper I, [Shettv. fc Ostrikeii 
I2006| ). and_jKim, & Ostrike^ II2002I ) for numerical examples, and 



are in fact natural consequences of "wiggle" instability, 
which is caused by a purely hydrodynamic phenomenon, 
i.e. Kelvin-Helmholtz instability (Wada, & Koda 2004|) 
(hereafter Paper I) . This phenomenon was also found in 
IShettv. fc Ostrikeii ()2006l ). They also pointed out that 
the features only grow in the inner-most several kpc re- 
gions if self-gravity and magnetic fields are ignored. 

More recently, three-dimensional (3D) response of the 
gas to the spiral potential was modeled using a local 
shearing box approximation in isother mal, MHD simu- 
lation s taking into account self-gravity (jKim fc Ostriked 
l2006l ). and found that the wiggle instability is suppressed 
by radial flapping motion of the shock. 

These previous results suggest that hydrodynamic ef- 
fects, self-gravity of the gas and magnetic fields play some 
important roles on the gas structures in a spiral poten- 
tial. However, an important feature of the real ISM has 
been ignored; its inhomogeneous multi-phase structures. 
Apparently the ISM is not 'isothermal fluid', neverthe- 
less it is assumed in most HD and MHD simulations'^. 
Effects of self-gravity and magnetic fields highly depend 
on gaseous temperatures and phases. In this sense, an 
energy equation with a realistic cooling and heating pro- 
cesses should be solved before taking into account those 
effects. Local box- approximation with a shearing peri- 
odic boundary (e.g. lKim fc Ostrikeii l2006f ) is not neces- 
sarily relevant for representing dynamics of the multi- 
phase ISM in galactic disks, because typical scales of 
the inhomogeneous structures of the IS M are not small 
enough compared to the disk size (cf. Wada. fc NormanI 
120011(20071 : [Tasker fc Brvanll2006[ ). Moreover, since grav- 
ity is a long-range force, global, 3D simulations for the 



ILa Vigne et al.i (i2006i) for observations. On the other hand 
'blanches' are sub-structures bifurcated from main spiral arms. 
They are smoother and longer azimuthal extent than spurs, 
which may be caused by resonances llChakrabarti, Laughlin. fc Shul 
[2003) . 



IDobbs. fc Bonnelll l|2007l ) recently studied gas dynamics in a 
spiral potential, taking into account multi-phase nature of the 
ISM. However, in their non-self-gravitating SPH simulations, an 
energy equation with realistic cooling and heating processes is not 
solved, alternatively warm (10"' K) and cold (100 K) components 
are treated separately as isothermal gases without phase exchange. 
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whole disk are essential if self-gravity of the gas is con- 
sidered. This is also necessary for the multi-phase ISM, 
because the scale height and velocity dispersion are dif- 
ferent for cold and hot gases. In 3D hydrodynamic sim- 
ulations of self-gravitating gas disks with the radiative 
cooling, we should take into account energy feedback 
from stars, otherwise the cold gas disk becomes unreal- 
istically thin and the gravitational instability is strongly 
affected. 

Here we show, for the first time, 3D evolution of the 
ISM in a galactic spiral potential, taking into account re- 
alistic radiative cooling and energy feedback from stars, 
especially from supernovae explosions. The simulations 
are performed for the whole gas disks without any as- 
sumptions for symmetry. In order to ensure a high spa- 
tial resolution (10 pc), which is essential to reproduce the 
multi-phase ISM, we here focus on a central part of a rel- 
atively small galaxy (radius is 2.56 kpc and the maximum 
circular velocity is 150 km s~^). However, this is large 
enough to see the global effect of the spiral potential on 
the inhomogeneous ISM. Our simulations clarify appar- 
ent discrepancy between the steady solution of "galactic 
shock" and the complicated, non-steady structures of the 
ISM in real spiral galaxies. This is a major progress from 
the previous simulations, and it will be a starting point 
for more realistic numerical models of the ISM, taking 
into account 'lived' stellar potential, magnetic fields, UV 
radiation from massive stars, and chemistry of molecules 
and atoms. 



2. METHODS 

Evolution of a gas disk (total mass is 4.3 x IO^Mq, ra- 
dius is 2.5 kpc) in a fixed gravitational potential <I>oxt is 
solved by an Eulcria n hydrod ynamic method with a uni- 
form Cartesian grid (IWada. fc Norman 2001; .Wadall2001l : 
IWada. fc Normanll2007f ). Here we brieflv summarize the 
numerical scheme. We solve the following conserva- 
tion equations and Poisson equation in three dimensions 
using 512 X 512 x 64 grid points with 10 pc resolution. 



| + v.(H = o, 



(1) 



— + {v-V)v + ^ + V$cxt + V*,g = 0, (2) 
ot p 



djpE) 
dt 

V2$,g = 47rGp, 



-V ■[{pE+p)v]^ pTuy + T,^ p'A{Tg), (3) 



(4) 



where, p,p, v are the gas density, pressure, velocity of the 
gas. The specific total energy E = |f P/2 + p/{'y — l)p, 
with 7 — 5/3. 

A cooling function A{Tg) (10 < Tg < lO^K) with 
solar metallicity, photo-electric heating by dust (Fuv) 
and a uniform UV radiation field (Go = 1.0) are as- 
sumed. The hydrodynamic sc heme is AUSM (Advection 
Upstream Sphtting Method) (|Liou. M.. fc SteffenI [l993h 
with MUSCL (Monotone Upstream-centered Schemes for 
Conservation Laws). The Poisson equation is solved to 
calculate self-gravity of the gas using the Fast Fourier 
Transform (FFT). 

The external potential is given by $cxt(^, 0, ^) = + 



$1 -I- $2: 

$o(i?, z) EEai;2(27/4)i/2(^2 ^ ^2 ^ z^^^^, (5) 
$i(i?, z) =ai;2 (27/4) 1/2 (i^2 ^ ^2 ^ z^)-^/^, (6) 
$2(i?, ^, z) = eobR^^i/iR^ + b^ + z2)3/2 

xcos[2(/) + 2coti •ln(i?/i?o)], (7) 

where i is the pitch angle {i — 10°) of the spiral potential, 
and i?o is an arbitrary constant, and a = 0.2 kpc, b = 2.5 
kpc, Va = Vb ~ 150 km s^^, and eo = 0.1. The pattern 
speed of the spiral potential is assumed Qp — 30 km 
s~^ kpc""'^. For comparison, we also run models with 
rip — 15 and 60 km kpc~^ and Eq = 0.05. However, 
we find no essential differences on the conclusions (see 
also discussion in §4). 

We take into account two feedback effects of massive 
stars on the gas dynamics namely: stellar winds and su- 
pernova explosions (F*). We first identify cells which 
satisfy criteria for star formation. The criteria for each 
cell where star formation is allowed are: 1) The gas den- 
sity is greater than a threshold, i.e. {Pg)i,j > Pc', 2) 
The temperature is less than the critical temperature, 
i.e. {Tg)ij < Tc, and 3) The criteria 1) and 2) are 
satisfied for 10^ yr. Here we have chosen Tc — 100 K 
and Pc = 600 cm^'^. Assuming the Salpeter Initial Mass 
Function with m„ = I2OM0 and mi — 0.2Mq, we cre- 
ate test particles representing massive stars (> 8Mq). 
Typically a few massive stars are replaced by one test 
particle. The initial velocity of each test particles is 
taken to be the same as its parent gas. The kinemat- 
ics of the test particles in the external potential and the 
self-gravitational potential of the gas are followed by the 
second-order time-integration method. The stars (test 
particles) inject energy due to stellar winds to the cells 
where the stars are locate d during their lifetim e which 
is approximately^ 10'' yr (Leith erer et al.iri992t l. When 
the star explodes as a supernova an energy of 10^^ ergs 
is injected into the cell where the test particle is lo- 
cated. Non-spherical 3D propagation of blast waves in 
the multi-phase, non-uniform ISM is then followed con- 
sistently in the hydrodynamic simulations. 

3. RESULTS 

Figure [T] is time evolution of the density field on x — y 
and X ^ z planes of the disk in the spiral potential. High 
density gases form two-arm spirals, and there are also 
substructures, i.e. "sp urs", between th e spirals, as seen 
in real spiral galaxies (|Elmegreenl 1 19801 ) . Each spur is a 
non-steady structure, but spurs are always seen in our 
simulations (See also the Supplemental Movie). As dis- 
cussed in previous 2D simulations (Paper I) each spur 
originates in a clump, whose size corresponds to observed 
a Giant Molecular Cloud, formed in the spiral arm. A 
small gradient of angular momentum in the clump is 
eventually enhanced, and the clump is stretched into the 
inter-arm region due to galactic rotation (see also Fig. 
[3|). In the present model, the wiggle instability (Paper 
I) is coupled with gravitational and thermal instabili- 
ties. This causes complicated distribution of the mat- 
ter coupled with the galactic rotation and tidal shear 
motions. In spite of these complexities, pseudo-spiral 
patterns roughly associated with the spiral potential are 
always present. Energy feedback from SNe (average su- 
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pernova rate is ~ 0.2 yr^^) change local density and tem- 
perature of the ISM, but it does not significantly affect 
global structures of the spirals in the galactic plane. Note 
that despite the complexity, the statistical structure of 
the density field (probability distribution function of den- 
sity) is represented by a single log-no rmal function over 
a wid e density range as predicted bv lWada. fc NormanI 

(Iloo3). 

On the other hand, as shown in the right panels of Figs. 
[T]and[21 vertical distributions of density and temperature 
are strongly influenced by SNe. Cold gases {Tg < 100 K) 
as well as hot gases (Tg ^ 10^^^ K) are blown out from 
the disk plane to the halo, where warm (Tg ~ 10''"^ 
K) media occupy most of the volume. The cold, dense 
gases sometimes form loop-like structures and fall back 
to the disk. In the inner region r < 1 kpc, where a su- 
pernova rate is higher than the outer region, the halo gas 
is occupied mostly by the hot gases. The interaction be- 
tween the disk and halo revealed here is a numerical rep- 
resent ation of a so-called "galactic fountain" or "chim- 
ney" llShapiro. fc Fieldl [1979 : iNorman. fc Ikeuchil 119891) . 
In the original idea of the galactic fountain, hot gas is 
spouted up into the halo from the disk, and it cools and 
subsequently falls back to the disk. However, our results 
show that not only hot gases, but also cold, dense gases 
are spouted up. The vertical structure is not steady, 
but in a long time-scale, the cold dense gas is distributed 
within ~ 100 pc from the galactic plane. The scale height 
of the gas is smaller in the central region due to the 
deeper gravitational potential, but effects of a spiral po- 
tential to the vertical structure are not clearly seen. 

In Fig. [21 a velocity field in a disk plane {z — 0) 
at t = 206 Myr is shown by streamlines overlaid on an 
iso-density surface^. The streamlines are bent near the 
spiral arms, and they are spread into the inter- arm re- 
gions. One should note that the "spurs" at the inter-arm 
regions are not located along the streamlines, but rather 
are perpendicular to them. As already reported in 2D 
simulations (Paper I), the spurs move along the galac- 
tic rotation, but they are stretched by the spread flow. 
This implies that the observed "spurs" are not waves. It 
is also apparent from Fig. [3|that the density and kine- 
matic structures are highly irregular, therefore it may not 
be represented by a local computational box with a pe- 
riodic boundary. Global simulations, in which the whole 
galactic disk is calculated without assuming symmetry, 
are essential to understand the structure and evolution 
of the ISM on a kpc scale. 

The spiral structure seen in Figs. [1]- [3] does not look 
like the typical hydrodynamic shock, in a sense that there 
is no clear jump of density, and it is consisted of non- 
steady substructures. Figure [4| demonstrates long-term 
evolution of density at r = 1.2 kpc. We see that hydro- 
dynamic shocks are formed upstream (i.e. on the concave 
side of the spiral) of the minimum of the potential trough 
in the first 20 Myrs or so. This is indeed the galactic 
shock solution found in the 60s for a tightly wound spiral 
potential. However, these shocks do not stay at the initial 
positions, and move back and forth between downstream 
and upstream of the potential minimum until t ~ 60 
Myr. After this early oscillating phase, the "shocks" dis- 

^ The isodensity-map is a method to represent a 3D density field 
by a opaque surface on which the volume density is the same. 



appear, and a new state emerges, characterized by its 
stochastic nature. The standing pattern near the poten- 
tial minimum is pseudo-spiral consisting of many clumps 
or parts of filaments. Each inclined narrow pattern seen 
in Fig. [4| represents a dense region orbiting circularly 
around the galactic center. The trajectory can be fol- 
lowed between the two spirals, i.e. for about 1/4 ^ 1/2 
of the galactic rotation. This means that the substruc- 
tures, i.e. spurs, in the inter-arm regions are linked with 
the non-steady 'chaotic arms'. Contrary to the classic 
galactic shock solution, the chaotic arms found here are 
located downstream of the potential trough on average, 
suggesting that the ISM does not behave like a smooth 
fluid. Typically the width of the waves is At/) ~ tt/6. 
Apparently the quasi-stable structures found here are no 
longer 'shocks' in a hydrodynamic sense. 

Finally, positions of massive "stars" overlaid on a den- 
sity map at t = 155 Myr are shown in Fig. 6. It is 
clear that massive stars form clusters, most of which 
are distributed near the potential trough, mostly down- 
stream, and they form spiral features. This is reason- 
able since the stars originate in cold, dense clouds. In 
real galaxies, the massive star complexes should form 
H n regions around them, and they illuminate the spi- 
ral structures. Our results show that the H n regions 
are offset to the downstream of the dust spirals, which 
is consistent with obse rvations (jdel Rio, fc Cepal Il998t 
iGittins. fc CTaTy [200l . However, as seen in Fig. 5, 
the spiral arms are not uniform, and positions of den- 
sity peaks in the pseudo-arms are not uniquely defined. 
Therefore, it would be difficult to derive dynamical infor- 
mation such as the pattern speed of the spiral potential 
and the position of corotation from observational data, 
for example the offset between B -band and Tband arms 
or the dust -arms and Hq arms (Idel Rio, fc Cepal [19981 : 
IGittins. fc Clark& .2004) . 

4. DISCUSSIONS 

4.1. Effect of the Pattern Speed of the Potential 

The density structure of the ISM in spiral po- 
tential can be affected by two free-parameters here: 
the pattern speed and strength of the spiral poten- 
tial. As mentioned in §1, bifurcation features of 
spiral arms could be explained by Lindblad and ul- 
traharmonic resonances (Shu. Milion e. fc Roberts 1119731 : 
IChakrabarti. Laughlin. fc Shu .2003 ). This substructure 
of the spir al shocks (i.e. "branches" ) is also seen in some 
models in lShettv. fc Ostrikeil (|2006l ). but we do not see 
clear evidence of the resonant structures in our results 
regardless of presence of resonances. Fig. [BJ; is density 
distribution at t = 61 Myr of a model with the pattern 
speed Hp — 60 km s^^ kpc~^ ^. As same as the fiducial 
model (Fig. [BId), in which flp = 30 km s~^ kpc~^, the 
density field is characterized with inhomogeneous spiral 
arms and spurs in the inter-arm regions. This is also the 
case in a model with flp = 15 km s~^ kpc^^ (there are no 

^ Since our disk is small and compact (the core radius Rcore is 
0.2 kpc and maximum circular velocity is 150 km s~^), Q{Rcore) — 
K/2(Rcore) ~ fl{Rcorti) /"i ~ Vmax / Rcor e is much larger than that 
in the l arge disk models assumed in lCha krabarti. Laughlin. & Sh^ 
(HoqD; IShettv. fc Ostrikeil iWOEi ). Therefore, the pattern speed 
assumed here is larger than th ose assumed in these previ ous papers 
(8.4 and 42 km s'^ k pc~^ inlShettv. fc Ostrikeii l|2006l), 21 .5 and 
11.5 km s^i kpc-i in lChakrabarti. Laughlin. fc Shul ^UM))- 
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resonances in the disk except an ILR at i? = 0.05 kpc). 
In Fig. [6K, the strength of the bar is half (i.e. £o — 0.05) 
of the fiducial model (Fig. (Hb) with the same pattern 
speed. Same as the fiducial models, the complicated spi- 
ral arms arc formed, but it is located more downstream 
than in the model with a stronger spiral potential at this 
moment. 

The average location of the spiral arms is probably af- 
fected by relative velocity of the gas to the spiral poten- 
tial and the depth of the potential. For smaller relative 
velocity (i.e. large Qp) and/or a deeper potential, the gas 
clouds tend to be trapped near the potential minimum. 
On the other hand, for smaller pattern speeds and/or a 
shallower potential, the clouds passes the potential min- 
imum and decelerate. As a result, density peaks tend to 
be formed more downstream on average. 

4.2. For Further Reality 

In the present model, we assume a time-independent 
two-arm spiral potential with a constant pitch angle, 
whereas the resultant morphology of the gas is quite com- 
plicated and time-dependent. These features are qualita- 
tively consistent with near-infrared and optical observa- 
tions of various types of spiral galaxies (Seigar,__& James 
Il998t iGrosbol. fc PatsisI Il998f) . In spite of the floccu- 
lent appearance of dust and young stars, the old pop- 
ulations seen in the K-hand are dominated by two-arm 
spirals. The pitch angles of K-hand arms are distributed 
between 5 and 10° independent of the Hubble type 
(Seigar. fc James! 11998'). The phase-offset between the 
potential, young stars and dust-arms in the real galax- 
ies could be used for studying the physical origin of the 
spiral arms. Unfortunately, observational constraints on 
the offset of the spiral arms of diff erent components are 
ambiguous (del R io, fc Cepal Il998f l. partly because the 
stellar and gaseous arms are not smooth, and it is hard 
to trace their positions, as is the case in our simulations. 
The distribution of the gas and massive stars are far from 
regular, but more complicated than the prediction of pre- 
vious idealized models. 

For further reality, one should take into account the 
non-linear gravitational coupling between the ISM and 
stellar potentials as well as the realistic treatment of the 
ISM itself. This should be achieved in self-consistent, 
cosmological simulations of formation of spiral galaxies. 
One shoul d realize, however that everi in recent simula- 
tions (e.g. ISpringel fc Hernguistj l2005l : iGovernato et al.l 
[2Q03), the mass and spatial resolutions for gas are 
1Q^~^Mq and sub-kpc, and therefore it is hard to repre- 
sent the complicated multi-phase structures of the ISM 
in spiral galaxies as seen in the present simulations. A 
much higher numerical resolution is also essential to re- 
solve the density waves in a stellar disk. 



It is not clear how the magnetic fields affect the sta- 
bility of the multi-phase ISM in a spiral potential. We 
would like to consider this problem in future 3D MHD 
simulations, in which realistic cooling/heating processes 
as well as self-gravity of the gas are at least considered. 

5. CONCLUSION 

The high resolution hydrodynamic simulations, taking 
into account self-gravity of the gas and realistic cool- 
ing and heating processes for the ISM, reveal for the 
first time that the galactic spiral arms of the ISM are 
neither hydrodynamic shock waves nor an assembly of 
long-lived, bullet-like cloudlets. The global spiral arms 
are consist of complicated time-dependent substructures 
from which stars can be formed, but over a long time 
(at least 5-6 rotational periods), they stably exist under 
the influence of the spiral potential. The pseudo-spiral in 
the multi-phase ISM is robust for energy input from the 
supernovae, which mainly cause the vertical non-uniform 
structure of cold and hot gases. The pattern speed of the 
spiral potential and its strength are not key parameters 
to alter above features. 

The ISM in galactic disks has often been represented 
by isothermal, non-self-gravitating fluid or inelastic par- 
ticles in many astrophysical simulations. Moreover, in- 
troducing a periodic boundary condition and reducing 
spatial dimensions were usual. However, those kinds of 
simplification do not necessarily represent the nature of 
spiral galaxies. Magnetic fields are not considered in the 
present simulations, but the present results suggest that 
even if the magnetic fields are weak, the spiral patterns 
can steadily exist on a global scale. Effect of magnetic 
fields and other important physical processes, such as UV 
radiation and chemistry, should be investigated based on 
3D global models, as presented in this paper. 

The present treatment of the multi-phase ISM could 
be used for direct comparison with observations, cou- 
pled with radiative tran sfer calculations for various 
observational p robes (e.g. IWada. Spaans. fc KimI 120001 : 
IWada. fc Tomis aka 200^. Fine structures of molecular 
gas associated with the spiral arms and in the inter-arm 
regions in external galaxies, which could be compared 
with the present numerical model, will be revealed by 
high resolution observations, e.g. by the Atacama Large 
Millimeter/Submillimeter Array. 
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Fig. 3. — Same as Fig. [T] but for streamlines on the galactic plane {z = 0) overlaid on an iso-density surface (p > 10 " '^^Mq pc ^ ~ 10 
cm""^) at t = 206 Myr. Note that the streamlines are generated by u and v components. The velocity field is non-steady, therefore the 
structure of the streamlines are time-dependent. 
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Fig. 4. — Time evolution of azimuthal density profile at r = 1.2 kpc. The color represents log-scaled density {Mq pc ^). Tfie white lines 
show the positions of the minimum of the spiral potential. The initial 'shock' turns to be stochastic arms consisting of many substructures 
after t ~ 60 Myr. 




Fig. 5. — Position of 7.9 x 10'* "massive star" particles overlaid on the density map at t = 155 Myr. The size of the box is 5.12X 5.12 
kpc. 
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Fig. 6. — Density distributions at t = 61 Myr for models with (a) Clp = 30 km s ^ kpc ^ and e = 0.05 (b) Hp = 30 km s ^ kpc ^ and 
e = 0.1, (c) Clp = 60 km kpc~^ and e = 0.1. Gray-scale represents log-scaled density in the unit of Mq pc~^. Radii of corotation, 
inner Lindblad resonances, and 4:1 ultraharmonic resonance are 2.7, 0.09, 0.74, and 1.46 kpc, respectively for (c), and 5.1, 0.05, 1.15 
kpc, respectively for (a) and (b). There are no ultraharmonic resonances in the computational region for (a) and (b). There are no 
clear resonance-driven features (so-called blanches) at any resonance radii. There is no clear sign that time-dependent inhomogeneous 
substructures of spiral arms and inter-arm regions (so-called spurs) are caused by the resononances. 



